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ABSTRACT 



Context. The Navarro-Frenk-White (NEW) density profile is often used to model gravitational lenses. Eor low values of the charac- 
teristic convergence (/c., « 1 ) of this model — corresponding to galaxy and galaxy group mass scales — a high numerical precision 
is required in order to accurately compute several quantities in the strong lensing regime. 

Aims. To obtain analytic solutions for several lensing quantities for circular NEW models and their elliptical (ENEW) and pseudo- 
elliptical (PNEW) extensions, on the typical scales where gravitational arcs are expected to be formed, in the Ks ^ I limit, establishing 
their domain of validity. 

Methods. We approximate the deflection angle of the circular NEW model, deriving analytic expressions for the convergence and 
shear for the PNEW and ENEW models. We obtain the constant distortion curves (including the tangential critical curve), which are 
used to define the domain of validity of the approximations, employing a figure-of-merit to compare with the the exact numerical 
solutions. We compute the deformation cross section as a further check of the validity of the approximations. 

Results. We derive analytic solutions for iso-convergence contours and constant distortion curves for the models considered here. We 
also obtain the deformation cross section, which is given in closed form for the circular NEW model and in terms of a one-dimensional 
integral for the elliptical ones. In addition, we provide a simple expression for the ellipticity of the iso-convergence contours of the 
pseudo-elliptical models and the connection of characteristic convergences among the PNEW and ENEW models. 
Conclusions. We conclude that the set of solutions derived here is generally accurate for Ks < 0. 1. Eor low ellipticities, values up to 
/f, - 0.18 are allowed. On the other hand, the mapping among PNEW and the ENEW models is valid up to /c, ^ 0.4. The solutions 
derived in this work can be used to speed up numerical codes and ensure their accuracy in the low /f, regime, including applications 
to arc statistics and other strong lensing observables. 

Key words, gravitational lensing; strong - galaxies: halos - galaxies: groups: general - dark matter 



1. Introduction: 

Gravitational arcs are powerful tools to probe the mass distri- 
bution in galaxies dKoopmans et al J [20091 : iBam abe et al. 201 1; 



Miralda-Escude 
2006h 



ft 



Suvu et al.l 120121) and clusters of galaxie s (Kovner 1989: 



20 lot iGilbank et al.1 1201 U I Wen et al.l l20Tll iMore et aP 1201 2t 
BavlissI 1201 21 : IWiesner et al.l 1201 2l Erb en et al., in prep), 
as well as in imag e s targeting clusters (iLuppino et al.l 19991: 
Eb eling et all I2OOII: IZaritskv & Gonzalez! 120031: ISmith et all 



T99I iHattori et al.l [T9971 IComerfordet aP 
In addition, their abund ance can he l p to constrain cos- 

I2OOI 



mological model s (Bartelman n et al.l 11998 : 



Bartelmann et all [20031: 
20101) 



iMeneghetti et al 



Golse et al.l 
I200I 



Julio et ail 



2005l:ISand et al.l2005i:lHennawi et alj20'08 l:'Kausch et alJ2010| : 
Horesh et al.ll201 U iFurlanetto et al.ll2012t i Postman etal.l|2012 ) 
and galaxie s (Ratnatunaaet al. 1999', Fassnacht et al. 2004| 



Bolton etalJ 12006: .Willis et al.. .2006: .Moustakas et al 2007 
Kubo & Deir Antonio! I2008t iFaure et al.ll2008t lJacksoi?l2008l) . 



Two techniques have been used to extract information from 
gravitational arcs. The first is arc-statistics: counting arcs as 
a function of their properties, such as the len gth-to-width ra- 
tio or angular separation , for a l ens sample dWu & Hamme 
1993!: iGrossman & Sahal 11994 iBarteknann & Weiss! 1199* 



gn 

3 



Bartelmann et al.! Il995t iBartelmannl Il995|) . The second is in- 
verse modeling: using arcs in individual clusters or galaxies aim- 
ing to dete rmine the mass d i stribution of the lens and source 
properties dKneib et al.lll993l:lKeetonll2001bLlGolse et alJl2002t 



IWavth & Webstei!l2006l:1jullo et al.ll2007ll2010!) 

These approaches have motivated arc searches in wide field 

I •— ' ■ 1 I H 1 I 

surveys^Gladders et al. 2003; Estrada et al. 2007; Cabanac et al 
12007!: iBelokui-ov et al.. .2009: .Kubo et al.i i2010: .Kneib et al 



Moreover, the upcoming wide-field surveys, such as the Dark 
Energy SurveyQ(Annis et al. 2005, Frieman et al, in prep; Lahav 
et al, in prep), which started taking data in 2012, are expected 
to detect of the order of thousands strong lensing systems, about 
an order of magnitude more than the current homogeneous sam- 
ples. 

A widely used model to represent the radial distribution of 
dark matter, from galaxy to clu ster of galaxies mass scales, is the 
Navarro-Frenk-White profile (iNavarro et al.lll996! 11997! here- 
after NEW) whose mass density is given by 



P(r) 



(r/r,)(l+r/r,)2' 



(1) 
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where Fs and p, are the scale radius and characteristic density, 
respectively. It is useful to define the characteristic convergence 



r^Ps 



(2) 



as a mass para meter, where, Scrit is the critical sur- 
face mass density ("Schnei der et al.l 119921: iPetters et all 1200 it 
[Mollerach & Roulet 2002). 

Some observational properties of many arcs systems (such 
as arc multiplicity, relative positions, morphology, etc.) imply 
that the mass distribution of the lens is not axially symmet- 
ric. Further, results from N-body simulations predict that dark 
matter halos are ty pically triaxial in s hape and can be m od- 
elled by ellipsoids dJing & Sutol 120021: iMaccio et al]l2007h . A 
first approximation to model realistic lenses is to consider el- 
liptical lens models, where the ellipticity is introduce d either on 
the m ass distribution (the so-called elliptical models. 'Sc hramm! 
[T99Q- Barkana 1998: Keeton 2001b: Ogurietal. 2003i) or on 
the lensing potential (the so-called pseudo-elliptical models, 
Blandford & K ochanek'19871: iKassiola & Kovnej[T99^ : iKneid 
20021: iGolse & Kneib 2002). While elliptical models are gener- 



ally more realistic, they are usually much more time consuming 
for lensing calculations than pseudo-elliptical ones. As a conse- 
quence, both have been used in the literature, depending on the 
application. 

In this work we consider elliptical as well as pseudo- 
elliptical lens models in which the radial mass distribution 
is given by the projected NFW profile. When used to rep 



resent lenses on galactic mass scales (see e.g., Asanol 

T^ — :^ ^4- .,1 lloAm'. '\7 4-*i p, T/- .^ lloAAnl. Ic.,,,., ^*- ^i I 



Davis et alj|20 03': Veg etti & Koopmansll200^ : ISuvu et alj 



200C; 



2013 



Ludlow et al.1 '2013) the characteristic convergence of the NFW 



model (Eq. (|2|) takes very small values. For instance, lensing 
systems with M200 < IO'^'Mq/i^^ Zl :^ 1 and sources with 
Zs = 2zl, have values o^ Ks < 0.05. In this regime a high 
numerical precision is required to accurately compute the func- 
tions involved in gravitational lensing, such as the deflection an- 
gle and its derivatives. Therefore, the numerical codes become 
either slower or, worse, provide unreliable results as we go to 
smaller values of k,. 

This issue has apparently not been addressed in the literature 
so far. A solution is to obtain analytical expressions, which are 
valid in this regime, providing at the same time a fast and reliable 
way to compute the relevant lensing functions. In this work, we 
present approximations for the deflection angle, convergence, 
and shear of the circular, pseudo-elliptical, and elliptical NFW 
models in the strong sensing regime. We use them to derive an- 
alytical solutions for iso-convergence contours, critical curves, 
and constant distortion curves. We compare these solutions with 
the exact calculations determining a domain of validity of these 
approximations. Further, for applications to arc statistics, we ap- 
ply these solutions to the calculation of the deformation cross 
section. 

The outline of this work is as follows: In Sect.[2]we present a 
few basic definitions of lensing quantities and introduce the no- 
tation used in this paper. In Sect.|3]we show the approximations 
for the lensing functions of the circular NFW model and their 
extension to the pseudo-elliptical NFW (PNFW) and elliptical 
NFW (ENFW) models. We also derive analytical expressions 
for iso-convergence contours and critical curves for these mod- 
els. In Sect. H] we obtain analytical solutions for constant distor- 
tion curves and for the deformation cross section. In Sect. |5]we 

^ To obtain this value of k, we use the expressions in Caminha et al. 
(in prep) with the same choices for NFW and cosmol ogical models as 
described in Sect. 2.2 of lDumet-Montova et al.l ( [2012f) . 



determine a domain of validity of these solutions in terms of the 
characteristic convergences and ellipticity parameters. In Sect.|6] 
we obtain a mapping among the PNFW and ENFW model pa- 
rameters. In Sect.|7]we present the summary and concluding re- 
marks. In AppendixlAlwe present the expressions for the poten- 
tial derivatives of elliptical models. In Appendix iB] we provide 
fitting functions giving upper limits of /c, for the validity of the 
analytical solutions as a function of ellipticity. 



2. Definitions and Notation 

We present below some basic definitions of gravitational lensing 
in order to set up the notation through out this work. More de - 
tails on the subject can be found, e.g., in'Schnei der et al.l (|1992|) : 
iPetters et al. (2001): Mollerach & Roulet (2002). 

The lensing properties are encoded in the lens equation, 
which relates the observed image position ^ to a given source 
position T] (both with respect to the optical axis) . By defining 
the length-scales ^0 on the lens plane, and 770 - ^oDqs/Di^ on the 
source plane, where Dqs and Dql are the angular-diameter dis- 
tances to the lens and source planes, respectively, the lens equa- 
tion is given in its dimensionless form by 



• 0!(X), 



(3) 



where _y - 1///70, x - ^/^o and a{x) - (^^^^^]a{^ox), where a 
is the deflection angle due to the lensing mass distribution. 

Properties of the local mapping are described by the Jacobian 
matrix of Eq. (O 



Jijix) = 6ij - diOjix). 

The two eigenvalues of this matrix are written as 



(4) 



Ai-(x) - I - k(x) + y(x) and At{x) - I - k(x) - y(x), (5) 
where 



k(x) = - (diaiix) + d2a2ix)) 



(6) 



is the convergence and y(x) - Jyl(x) + ylix), is the shear, 
which has components 

nix) = - (diaiix) - d2a2{x)) > 72 = r ^^^"'^{x) + d2ai(x)) .(7) 

Magnification in the radial and tangential directions are given 
by /1,t'(jc) and Ay^{x), respectively. Points satisfying the condi- 
tions Ar{x) - and A,{x) = are the radial and tangential crit- 
ical curves, respectively. Mapping these curves onto the source 
plane, we obtain the coiTesponding caustics. 



3. Solutions for iso-convergence contours and 
critical curves 

In this section we present the approximations for the lensing 
functions of the the circular NFW (Sect. ITTJ, the PNFW (Sect. 
O, and the ENFW (Sect.O models. 



3.1. Circular NFW model 
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3.2. Pseudo-Elliptical NFW model 



Following iBartelmannI (Il996l) . from the density profile dD and The construction of pseudo-elliptical models is made by replac- 
ing the radial coordinate of the lensing potential by 



taking ^o - fs, the dimensionless deflection angle is given by 



a{x) — 

X 



1 ^ 

2 vr 



: arctanh 



1-x 

1 +x 



(8) 



Xip = X Ja^p cos^ (p + btp sin" 0, 



(17) 



from which the convergence and shear are derived. 

In the limit of /c, <k 1 the typical scales corresponding to 



where a^ and b^ are two parameters defining the lens- 
ing potential ellipticity. Adopting the approach of the A ngle 
Defl ection method introduced by iGolse & KneibI (l2002h (see 



the strong lensing regime (e.g., the size of critical curve and also jDumet-Montoya et al. | HqH hereafter DCM), from Eqs. 

(I9b- (llll l. the deflection angle, convergence, and components of 
the shear of the PNFW model are 



caustics) are much smaller than the unity. In this case a high 
numerical precision is required to accurately compute lensing 



quantities, such as the shear and convergence. However, simple / i — . \ 

analytic expressions can be found in this regime avoiding such 0!ip(^) - o:(Xtp) ^ ya^ cos 0^, yfc^ sm (ptpj , 



numerical difficulties. Keeping the first terms in a series expan- 
sion of ([8]l for X <& I leads to 



a(x) - -xKs 1 1 + 2 In - 1 . 

In this case, the convergence and shear are given by 

'^^^-2[-x^i;r-H'^'^2y 

y{x) 



2\~ 



da 
dx 



(9) 

(10) 
(11) 



K^{x) = JIk(x^) - Sk"^ COS 2(f)^, 

yi^ix) = !BK(xy;) - JIk^ cos 20^, 

r2<p(-«^) = - Ja^b^K^ sm24>^, 
where 



(18) 

(19) 
(20) 

(21) 



Allowing a relative deviation of less than 5% with respect to the 
exact expression, we found that Eq. (|9]l is a good approximation 
to the deflection angle for x < 0.25, while the same holds for 
Eqs. ( [Tol l and ( fTTT i for the convergence and shear for x < 0.15. 

From Eq. (fTOt and fixing a value for the iso-convergence 
contour as /Cconst, the equation k{x) = /Cconst has the solution 



^ = :z(a^+b:p), S = -(a^-b^) and (p^ = arctan( J/j^/a^tan^j, (22) 

and we denote the NFW characteristic convergence, Eq. (|2}, by 

'fi 

From Eq. ( fT9] l and fixing a value for the iso-convergence 
contour as /Cconst, the equation k^{x) = /Cconst has the solution 



Xk(<P) 



y/jl + S COS 20 



/fconst + (2y[ + S COS 20^)/f^ \ 

expl ^^ |.(23) 



2exp 



'^const "'" ^^s 

2a: t 



Therefore, the iso-convergence contours are not elliptical, as is 
well known for pseudo-elliptical models. 

The solutions for critical curves are a bit more involving. 
(12) From Eqs. (fT9ll- (l2n i the determinant of the Jacobian matrix is 



From Eqs. ©, ([TOll and ^TT^ it follows that 



A,ix)^l+KJl+2\n^], Arix)=l+KJ3+2\n'^] (13) 



detJ(jc) = a^bipK^(Xip)-2JlK{x^)+l+2!BKt cos2(f)^-a^by;(K'^)^ .(24) 
Then, solving the equation det J(a:) = for k{x^), defining 



2/' '^ ' 'V 2) 

such that the determinant of the Jacobian matrix is 

det J(x) = K-{x) - 2k(x) +1-k1. 



(14) 



K^(4>) = (a^b^)-^ (jl ± ^& + (a^b^Ky- - 2a^b^SKt cos 20^) , 
and inverting K(x^p) - ^^(<^), using Eq. ( fTTl i. we obtain 



The solutions for the tangential and radial critical curves follow 
from the equations above. The radial coordinates of these curves 
ar^l 



x,(4>) = 

Xri(f>) = 



X, - 2exp(- 
Xr = 2exp 



1 + Ki 

2ks 

1 + 3A-.5 

2k ^ 



2 
^[W+Bcos24> 



: exp - 



exp 



^-(0) + 2/f? 

Tt{(l)) + 2£, 

2ll 



(25) 
(26) 



(15) 
(16) 



These curves are mapped onto the source plane by using the lens 
equation with the deflection angle given in Eq. ( fTSl ). 



3.3. Elliptical NFW model 

The expressions for the caustics are straightforwardly obtained ,,, , „-_„-, , ,, , . , ,. , ,. 

,_., : „,:„_ ,, : , :„ ^,. _ T ,: :,,. ,,. _ Wc construct the ENFW model by replacmg the radial coordi- 



nate of the surface mass density, Eq. ( fTOl i. by 



by inserting the expressions above in the lens equation with the 
deflection angle given in Eq. ^ 

3 These solutions are shown in iMeneghetti et al.l ( l2003h . However, x^ ^ x yjcos^- <pla\ + sin^ (pib 
there is a typo in this reference as they appear in their Eq. (11) as the 
eigenvalues of the Jacobian matrix. where a^, and b^, define the ellipticity of the mass distribution 



(27) 
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The lensing functions of this model can be written as (see 
AppendixlAt 



aiiix) - a^bxx^J'iKix) - 2i^£ii(/>)jsin(f>, 

a.b.(Jo + Ji)k(x)] , 
1 



K^{X) = K{X^)^1-^1'P+Q- 



mix) 



la^b^(Jo-Ji)K(x)+Q-P]., 



r2i:(-«^) = -flsZj^/Cj sin207Cii 



(28) 
(29) 

(30) 

(31) 
(32) 



where we denote the characteristic convergence, Eq. (|2), by k^ 
and define 



^ = 1 + 2flx^2^ [KoW cos2 + Xo(0)] , 
(3=1+ la^b^i^ [KiW sin- + i;i(<^)j . 



(33) 
(34) 



4. 1 . Constant distortion curves 

A typical arc forming region in the lens plane is determined 
by the so-called constant distortion curves, corresponding to 
the \Ra\ = Rth contours. An often used value for Rth is 10. As 
the value of R± is decreased, the inner curve (corresponding to 
R^ - -Rth) gets closer to the center, while the outer enclosing 
curve (Ra - +/?th) reaches higher radii, where the analytic ap- 
proximations derived in section 3 are less accurate. Therefore by 
using a lower value of Rth to determine the limit of validity of 
these approximations we are assuring they are even more accu- 
rate in the arc formation region. For this sake, we adopt Rth = 5 
when we make the numerical comparisons to the exact solution 
throughout this paper 

From the approximations given in Sect. |3] it is possible to 
obtain analytical solutions for the radial coordinates of constant 
distortion curves. For the circular NFW model, from Eq. (ITjI . 
the equation R,\{x) = Rth has the solution 



For an iso-convergence contour value /Cconst, the solution for 

K-l(x) = A-const is 



xa -2 exp 



4 



:exp - 



cos^ (p/al^ + sin (p/bi 



^const "•" ^^s 

'2^ 



2KARth-l) 



(40) 



(35) 



which are ellipses, as expected. 

From the definitions in (|5]l and Eqs. (I30ll-(l32li. the determi- 
nant of the Jacobian matrix is 

det J(;c) = {a^b^fjoJiK\x)-a^b^(Jir+JoQ)K{x)+rQ-yl^.(36) 

Following a similar procedure to obtain the critical curves for 
the PNFW model, the solutions of the equation det J(jc) = are 



For the PNFW model, the calculation of the radial coordi- 
nates of the constant distortion curves is a bit more complicated. 
From Eqs. (fT9]l- (l2n i. solving the equation R^ix) = R^ for k{x^), 
we obtain 



k(x^) 



were 



I^Aip 



2[m^-&Qiy 



(41) 



Xt{4i) = 2exp - 



Xri^i) - 2 exp - 
where we define 



K-^((t>) + 2KY 

2^? J 
2^ 



(37) 



(38) 



with 



■'Ri^ - ^(^1^ - Ki^f + 4(J?[2 - &Ql)Qlyl^ 



~4(4» 



Ji'P + JoQ± ^(Ji'P-JoQ)^+4JiJoyl^ 



2a^b^oJi 



Rth + 1 
Gth - r, 

"th - t 

^1^ = (^ - Seth)(l + (^gth + S)k'^ cos 20^), 
^2^ = (J?[ + Seth)(l - (^gth - S)k'^ cos 20^), 

where (p^ is given in Eq. (l22l i. Inverting Eq. Wi\ . using Eq. (fTTl i. 
we obtain for any angular position 



The corresponding caustics are obtained by using the lens equa- 
tion with the deflection angle given in Eqs. ( l28b and ( 



4. Solutions for the deformation cross section 

Gravitational arcs are usually defined as images with length-to- 
width ratio, L/W, greater than a threshold Rth- For fast calcula- 
tions in arc statistics, it is useful to approx imate L/W to the ratio 
of the eigenvalues of the Jacobian matrix ( Wu & Hammeill993t 
iBartelmann & WeisJl99"4tlHamana & Futamaselll997h 



^Aip ■ 



V^ + S cos 24> 



exp 



'<,iip 



4{je - &QI)k\ 



fP\ 



4(^2 - S^Q? )4 



(42) 



Following the same procedure above, for the ENFW model, 
from Eqs. (I30ll-(l32li. we obtain for each angular position 



xxi. = 2 exp - 



Kay + 4Q\y.Q2y.k] 



W ^ l^^l 



(39) 



where R^ - Art At. 

This approximation holds for infinitesimal circular sources 
and bre aks down for arcs generated by the merger of multiple 
images jRozo et al.ll2008h or by large or non circular sources. 

In this section, using the approximation above, we derive 
analytical solutions for constant distortion curves for the NFW 
models (Sect. l4Tb . We thereafter employ these solutions to com- 
pute the arc cross section (Sect. |4j2]i. 



4<3i2(32E^? 
where we have defined 

yjCRiYQiY-'RizQiY)- + leQi^Qi^Qlrl^ 

with 

!Riz = P + Q+QthCP-Q), 

HiY = r + Q-QthCP-Q), 

QxY = a^bAJt,+Ji+Qth{Ja-Jx)l 

QiY = a^b^Jt, +Ji- QthiJo - Ji)l 



(43) 



Dumet-Montoya et al.: Analytic Solutions for NFW lens models 



where V, Q and J'„ are given in Eqs. (l33T l. ( l34l i and (IA.8b . re- 
spectively. 

In Eqs. ( I40b . ( l42b and ( |43] ) the solution for the equation 
R,i{x) - -Rth is obtained by replacing /?th — > -^th- Also, at 
the limits R^h -^ oo and Rth — > 0, these expressions yield the ra- 
dial coordinates of the tangential (Eqs. (flST l. (IZST i and (l37T i) and 
radial (Eqs. (fTSI l. (|26] | and ( l38T l) critical curves of the correspond- 
ing models. 

4.2. Arc cross section 

The arc cross section, ctk^^, is defined as the weighted area in the 
source plane, such that sources within it will be mapped into arcs 
with L/W > Rth- This cross section is usually computed using a 
large sample of arcs obtained from ray-tracing an even larger 
number of finite sources and is computationally demanding. An 
alternative for fast calculations is to use the approximation ( l39l l. 
In this case, the arc c ross section is calculated in the lens plane 
as (iFedeli et alj|2006l DCM) 



o"«,„ 



CTp , cr 



\D0LJ 



o"R,h; CTRih 



/ 

J\RA\>Ra. 



\detJ{x)\xdxd(l), 



(44) 



where the quantity o-r,,^ is known as (dimensionless) deformation 
cross section. 

For small values of the characteristic convergences, the de- 
terminant of the Jacobian (see Eqs. (fT4l i. (l24l i and ( |36] |) takes the 
form 

det J(a:) = Ak^{x) + Bk{x) + C, 

where A, B, and C are independent of the radial coordinates. It is 
possible to reduce the calculation of (l44l i to a one-dimensional 
integral. Inserting the expression above into Eq. (l44t and inte- 
grating over the radial coordinate, within the lower and upper 
limits given by the constant distortion curves (i.e., from -Rth to 
Rth), we obtain 



1 r-" 



c^«,h = :^ I ['^(x+) + S(x^) - 2Six,)] d(p. 
For the circular NFW model we have 

S(X) = X^ [1 - K, - K(x)f , 

such that Eq. ( |45] ) gives 



4-KKi 



0"«,h 



(Ri-^)' 



[{Rth + i)'4 + (Ra 



dvi 



(45) 



(46) 



(47) 



where x± are given in Eq. ( l40t for Rth and -Rth, respectively. 
This expression shows the exponential dependence of the cross 
section on k, (Caminha et al., in prep). For instance, varying 
Ks from 0.01 to 0.1, os changes approximately by 41 orders of 
magnitude. Further, o"r,^ oc R^^ for Rth » 1 as expected from 
the be havior of sources near to the caustics (iMollerach & Rouleg 
I2OO2I Caminha et al., in prep.). 
For the PNFW model we have 



S((p) - x^ [a^b^S^i^ - (a^ + b^)S\^ + >So^j , 
where 

Sq^ - \ +{a^- b^)K'^ cos 2(f>^. 



(48) 



In this case, this function must be evaluated at x, and x± given in 
Eqs. (l25T l and (l42l i. respectively. 
For the ENFW model we have 

Sicf,) = x^ [(a^b^f{JoJi)Si^ - a^b^iJiP + JoQ)Si^ + Sq^] , (49) 

with 

This function must be evaluated at x,, Eq. (ISTT l. and x-t Eq. ( l43T l. 

5. Domain of validity of the solutions 

In this section we quantify the deviation of the analytic solutions 
(Sects. [3]and|4]i with respect to their corresponding exact calcu- 
lations, seeking to determine their domains of validity in terms 
of the model parameters (Sect. lSTt . Then, aiming to test the do- 
main of validity of these solutions for computing other lensing 
quantities, we compare the deformation cross sections in Sect. 



©^ 10-^ 




0.25 



Fig. 1. Mean weighted squared radial fractional difference £)- 
for the constant distortion curve with Rth = 5 for the circular 
NFW lens model as a function of the characteristic convergence 



5. 1 . Limits for constant distortion curves and critical curves 

To compare the analytic solutions for critical curves and con- 
stant distortion curves to their exact calculations, we use as a 
figure-of-merit the mean weighted squared fractional diff'erence 
between the two curves (Dumet-Montoya et al. 2013) 
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Fig. 2. Mean weighted squared radial fractional difference £)^, for constant distortion curves with /?th = 5, as a function of the 
characteristic convergences, for some values of ellipticity parameters. Left Panel: PNFW model. Right Panel: ENFW model. 



where A^ is the number of points of the curves, 0, is their polar 
angle, w,- = 0,- - 0,_i is a weight (to account for a possible non- 
uniform distribution of points), XaS.^i) and x^^^i^t) are the radial 
coordinates of the curves curves obtained from the exact and 
approximated calculations, respectively. Notice that, the D~ in 
Eq. (ISOl l is independent both on the lens length scale and on the 
discretization. 

Choosing a maximum value for £)^, we can define upper val- 
ues for the characteristic convergences (for a given ellipticiy) 
such that the contour curves obtained with both the exact and 
approximated calculations will be close enough to each other. 
This maximum value is chosen by visually comparing the ap- 
proximated and exact solutions for critical curves and constant 
distortion curves for several values of D^ and combinations of 
the NFW model parameters. 

We compute £)- on the 7J^ - +Rib curves since their points 
are the farthest from the lens center Thus, imposing a limit on 
the NFW parameters to match the R^ - Rih curve, we will auto- 
matically match the curves enclosed by it. We have checked that 
for a broad range of the NFW model parameters and ellipticities 
(for both elliptical and pseudo-elliptical models) a good visual 
matching of the exact and approximated curves if obtained for 
a maximum value of TP' = 10"'* (for R± = 5). This value also 
ensures a good matching of the critical curves and we found that 
the corresponding curves in the source plane are well-matched 
too. We will thus fix this as the upper value of D^ throughout 
this work. 

In Fig. [T] we show D~ for the /?th = 5 curve as a function 
of Ks for the circular NFW model. Setting the above upper value 
for D^ leads to a maximum value of /c^ - 0.18, for which all 
curves |^^| = Rthi^ 5) as well as the radial critical curves are 
well matched. 

For the PNFW model we chose the convention 
(iBlandford & KochaneklllQSTtlGolse & Kneibll2002h 



= 1 - e, bip — I + e. 



(51) 



where the ellipticity parameter s is defined in the range < 
e < 1 . As shown in DCM, for small values of k^ we must have 



E > 0.5 to avoid dumbbell-shape mass distributions. Hence, in 
the forthcoming analyses, we adopt an upper value of e = 0.5. In 
the left panel of Fig. |2] we show D^ as a function of k^ for some 
values of s. For the higher value of s considered, D^ - lO"'* 
is reached for /c^ = 0.1, providing an upper limit for the char- 
acteristic convergence for the validity of the approximations. 
However, higher values of k^ are allowed as e decreases. For 
instance, when e — > 0, we find that /Cj — > 0.18 (in agreement 
with the result for the circular NFW model). 

For the ENFW model, we chose the parameterization 
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fls 



^/^^ 



and by 
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■fix, 



(52) 



where the elUpticity s^ is defined in the range < ei < 1 . In 
the right panel of Fig. |2] we show D^ as a function of k^ for 
some values of e^. The behaviour of £)- as well as the maximum 
values of ez for each /cf are qualitatively similar to that of the 
PNFW model. In this case, we find that a maximum value of 
K^ - 0.1 is allowed for s^ - 0.7. Again the maximum value 
KJ -^ 0.18 as £i ^ 0. 

In Appendix|B]we present fitting functions for the maximum 
values of s {s-£) as a function of /Cj (k^) for which the approxi- 
mation is accurate following the criteria of this section. 



5.2. Comparison between deformation cross sections 

We compare the exact calculations and the solutions derived in 
Sect. 14.21 in order to verify if the limits derived in Sect. 15. II also 
give a domain of validity for the deformation cross section. To 
quantify this comparison, we compute the relative difference 



Ad-R 



o"«,i, 



o"R,i, 



-o-R, 



ih.app 



o-R,h 



(53) 



where the subscripts ex and app refer to the exact and approxi- 
mated calculations, respectively. 

In Fig. [3] we show 0-5 (upper panel) and Ao-^/ds (bottom 
panel) as a function of /c, for both the exact and approximated 
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Fig. 4. Deformation cross sections for the PNFW model (left panel) and ENFW model (right panel). Symbols correspond to the 
exact calculation. 
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Fig. 3. Deformation cross section for the circular NFW lens 
model as a function of /c,. Solid line corresponds to expression 
Wh . Filled circles correspond to the exact (numerical) calcula- 
tion. 



calculations. Considering the limit of the previous section, i.e., 
Ks < 0.18, Eq. (l47b deviates from the exact calculation by at most 
2.5%. 

In the left panel of Fig. |4] we show 0-5 (upper panel) and 
Aa^/as (bottom panel) as a function of /Cj for the PNFW 
model for both the exact and the approximated calculations. 
Considering the limit of applicability obtained in the previous 



similar result is found for the ENFW model. In this case, for 
0.1 < /fj < 0.18 and sy, < 0.7, we find that Ads/ 0-5 (o"5 calcu- 
lated by Eq. (|45| with Eq. (|49])) is at most 2.5%. 



6. Mapping among PNFW and ENFW parameters 

As mentioned in the introduction, elliptical models are more 
physically motivated than pseudo-elliptical ones. However the 
determination of lensing quantities, suc h as the s hear, for t he for- 
mer re quire the evaluation of integrals (ISchramm 1990: KeetonI 
l2001al see also Appendix [All, which generally have to be com- 
puted numerically. On the other hand, pseudo-elliptical mod- 
els do not require such integrals to be evaluated (see, e.g., Eqs. 
(fT8ll-(l2T])) allowing for fast calculations for the same quantities. 
However, it is well known that pseudo-elliptical models have two 
main limitations: their s urface mass density can assum e negative 
values in some regions (iBlandford & Kochanek"1987) and may 
present a "dumbbell" shape for high ellip ticities (Kovner 1983; 
[Schneider et"anil992l: iKassiola & Kovneilll993l) . At least in the 
case of the NFW model, the first problem does not affect the 
region of arc formation (DCM). Regarding the shape of the iso- 
convergence contours, for each value of k^ there is a range in s 
for which the contours are approximately elliptical (DCM). In 
principle, within this range of parameters, the pseudo-elliptical 
models could be employed instead elliptical ones in studies that 
require numerous evaluations of lensing quantities. For this sake, 
a correspondence among model parameters has to be estab- 
lished, associating a pair of the PNFW parameters (e, k^) to a 
pair of the ENFW ones [s^, k^). 

To obtain s^ from the PNFW model we use the same proce- 
dure as in Golse & Kneib (2002, see also DCM). From Eq. ^ 
we define the semi-major axis a by Xiti(p - 0) and, the semi- 
minor axis b by Xit{4> = n/2) such that 



section (0.1 < < < 0.18 and e < 0.5), Eq. 



with Eq. 



( |48] ), deviates at most 2.6% from the exact calculation. A very 
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= 1- 



1 -s 



1 



■ exp 



(54) 



Notice that this expression has no dependence on 7?th nor on /Cj . 

The quaUtative behaviour of ej;(e) from this equation is very 
similar to what was found numerically for generic values of /f^ 
in DCM (see e.g., their Fig. 6). For e «: 1, the expression above 
gives ££ ^ 2e, in agreement with DCM (see e.g., Eq. (B.2) which 
gives sx = 1.97e, for /c^ — » 0). 

Interestingly Eq. ( l54l i is a good approximation to the exact 
relation for values of /c-^f well above the limit of validity of the ap- 
proximation derived in Sec.|5] Indeed, the relative deviation with 
respect to the exact calculation of ej;(e, a-^, 7?th) (see the Elliptical 
Fit Method in Sect. 4 of DCM) is at most 5% for kI, < 0.4, with 
E < 0.5, and/?th > 5. 

To associate a value of k^ to a pair (a-^, e) we impose that 
the tangential critical curves of the PNFW and ENFW models 
match at = 0. This condition, from Eqs. dZSl l and (|37] |. gives 



/^ 






a^b^ [Ji + b^K'^iJi log a^-Ji- 2Xi (0))] ' 



(55) 



Note that, for e = e^ = the expression above reduces to k^ - 
K^ = Ks, as it should be. 

We determine the upper value of k^ such that the expression 
above deviates by at most 5% with respect to the exact calcula- 
tion. This yields k^ < 0.5 for s < 0.5 and Rth > 5. Again, the 
domain of validity of the approximate mapping is larger than for 
the lensing equations. 



7. Summary and concluding remarks 

When considering small values of the characteristic convergence 
of the NFW model (i.e., at the galactic and galaxy group mass 
scales), some strong lensing quantities require a high numeri- 
cal precision to yield accurate results, which can be time de- 
manding. Motivated by this issue, we obtained analytic solutions 
for several strong lensing quantities for elliptical and pseudo- 
elliptical NFW models and quantified their corresponding limits 
of validity. 

The starting point is the approximation (jS) for the deflec- 
tion angle of the circular NFW model. This approximation was 
applied to the standard prescriptions for obtaining the conver- 
gence and shear for circular, elliptical, and pseudo-ellitucal mod- 
els, leading to analytic solutions for these quantities (Sect. |3]l. 
Those were in turn used to derive analytic expressions for iso- 
convergence contours and critical curves (see Sect. |3) and for the 
constant distortion curves as function of R± (Sect. 14.11 ). 

As a practical application of these results, we computed the 
deformation cross section (ctr^i ). In the case of the circular NFW 
model, we obtained an analytical formula for ctr^^^ (Eq. (|47])). 
which reproduces the behavior with respect to k,, and the scal- 
ing with Rth obtained numerically in previous works (see e.g., 
Caminha et al., in prep). We have shown that the computation 
of this cross section is reduced to a one-dimensional integral for 
both the PNFW and ENFW models (Eq. ^ with either Eq. (|48]l 
or Eq. ( |49] l, respectively). 

We used the figure-of-merit D~, Eq. ( l50l ), to quantify the de- 
viation of the solutions of the constant distortion curves with 
respect to their exact calculations. Setting a maximum value of 
1)2 _ J 0^4^ we find that Eq. (gSJi (Eq. (|^), matches its corre- 
sponding exact calculation up to k^ = 0.1 {k^ = 0.1) for s < 0.5 



(fis < 0.7) and 7?th > 5. In particular, we find that the char- 
acteristic convergence can go up to 0.18 as the ellipticity pa- 
rameters tend to zero, as expected from the limit derived for the 
circular NFW model. In Appendix iBl we provide fitting func- 
tions for the maximum ellipticity allowed as a function of the 
characteristic convergences (in the range 0.1 < k^,k'^ < 0.18) 
to ensure the validity of the approximations. We emphasize that 
these limits, also ensure a good matching for critical curves and 
the Ra - -Rth curves to their corresponding exact calculations, 
since these curves are enclosed by the R^ - Rth curve. We veri- 
fied that the corresponding curves in the source plane also match 
the exact solution for the parameters within the limits derived 
above. 

We compared the deformation cross sections (obtained from 
both the exact and approximated calculations) in order to check 
if the domain of validity derived for the constant distortion 
curves also holds for this quantity (Sect. 15.2) . which we found 
to be the case. For instance, for the circular NFW model, for 
Ks < 0.18 and Rth = 5, Eq. ( l47T i deviates of at most 2.5% from 
the exact calculation, while for the PNFW and ENFW, Ao-s/cts 
(Eq. (|53]l) is at most 2.5%, for 0.1 < /c^ < 0.18 and e < 0.5, and 
0.1 < /f^ < 0.18 and s^ < 0.7, respectively. 

Overall, the approximate solutions presented here are accu- 
rate for all strong lensing quantities that were addressed in this 
work, within the considered parameter ranges, for /c, < 0.1. In 
some cases they are valid for higher /Cj, such as for low ellitpci- 
ties. Furthermore, some derived quantities are valid up to much 
higher characteristic convergences. For example, the ellipticity 
of the iso-convergence contours of the PNFW model is well re- 
produced by the simple analytic form of Eq. (l54t up to k^ ^ 0.4. 
We found that in this range the relation s^{s) is independent of 
K^ and the chosen value of the contour To complete the asso- 
ciation of the PNFW to ENFW model parameters we derived a 
relation among characteristic convergences, Eq. (l55t . which also 
matches the values obtained in DCM for /Cj < 0.5 and Rth > 5. 

The analytic solutions presented here allow for a robust and 
fast computation of several strong lensing quantities to be car- 
ried out in the low characteristic convergence regime. They may 
thus be useful for the lensing community and could be readily 
included in Strong Lensing codes, considering the domain of ap- 
plicability derived in this work. 

In principle, the approximate solutions derived in this paper 
could be extended to other quantities, such as the lensing mag- 
nification and the magnification cross section. They might also 
be applied to find solutions of the lens equation, for multiple 
images and arcs. One way of finding approximate solutions for 
arcs, for low ellitpicities, is through th e Perturbative Approach 
(lAlardll2007l 120081: iPeirani et al.ll2008l) . for which the Einstein 
ring solution (given analytically in Eq. (flST l) is the starting point. 
Therefore, analytic solutions for arcs might be derived in this 
framework. The investigation of these and other possible exten- 
sions is left for future work. 
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Appendix A: Lensing Functions for elliptical lens 
models 

The lensing functions of models with elliptical mass distribution 
can be obta ined following th e expressions derived in iSchrarnml 
d 19901) and iKeetonI ( 1200 lal) . which were generalized for any 
choice of the ellipticity parameterization in Caminha et al. (in 
prep.) and are given by 



an = aj^b-^Jox cos (p, 021: - ai,bj, J ix sine 
diaij, - aj,b-i:yJo + 2Kox^cos^ 4>j, 

820121. - ajjjj, (/i + 2K2X^ cos^ 4>) J 

dia2i. — a-£bY,K\xsm2<p, 

where aj; and b^. are defined in Eq. (IZTT i. 



J n •■ 

and 
K, 



K{m{u)) 



[1 - (1 - /7|)m]1/2+"[1 - (1 - 4)m]V2-« 



Jo 



UK'{m{u)) 



[1-(1 -bl)uYI^+"[l-(l -a2)M]5/2-„' 



(A.l) 
(A.2) 

(A.3) 
(A.4) 

(A.5) 
(A.6) 

1 dK{m) 



where k is the convergence of the circular model, k' - 2,,, j,„ 
and we have defined the variable trp- = x^g{u, (p) such that 

2 ( cos^ (p sin^ 

g (u, (f>) - u\ — + — 

U-(1-4)m 1 -(1 -bl)uj 

For the ENFW model, with the approximation (fTOl i. we can 
rewrite the potential derivatives as 



Jn(x, <f>) = J„K{x) - 2k';Z„{cP\ K„(x, 4,) = -^7C„(0), (A.7) 

X 



where 

jn = r 

Jo 

Jo 
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■Cnicp) = 



du 
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-/7^)m]1/2+«[1 -(1 - 

ug{<p, uy^du 


- al)u]y2-" 
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-(1- 


log g((p,u)du 


- 4)m]5/2-" 



(A.8) 
(A.9) 
(A. 10) 



[1 - (1 - /72)m]1/2+«[1 - (1 - 4)m]3/2-«' 

Substituting the expressions above in Eq. dA.lb we obtain 



Eqs. (|28ll and ^. Similarly, by substituting Eqs. (lAJJ i - ( lAlOl i 
in ( IA.2b - (I A. 4b and using the definitions ^ and (|7]) we obtain 
Eqs. (I30| - 01. 

Appendix B: Fitting functions 

Applying the procedure outlined in Sect. 15.11 we obtained the 
maximum values of e (ex) as a function of /Cj (/<^) such that the 
analytic solutions derived in the paper are accurate. We find that 
these functions are well fitted by a Pade approximant of the form 



e{Ks) 



2jn=0 '^n(Ks) 



(B.l) 



where e and /c, correspond either to e, a-^ (for the PNFW) or 
gj;,*^ (for the ENFW), and the values of coefficients are given 
in Tableim] 





PNFW 


ENFW 


ai 
as 


0.187 
-0.770 
-0.556 
-3.129 


0.300 
-1.236 
-1.481 
-0.796 


bo 
bi 
b2 


-0.024 
1.670 
5.021 


0.125 
0.542 
4.717 


X^ 


4.64 X 10-^ 


8.69 X 10-" 



Table B.l. Results from the regression analysis using the Fade 
approximant for ^(a:,). The last row corresponds to the values of 
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